arXiv:1501.07611vl [astro-ph.SR] 29 Jan 2015 


Draft version February 2, 2015 

Preprint typeset using DTeX style emulateapj v. 5/2/11 


PREDICTIONS FOR OBSERVING PROTOSTELLAR OUTFLOWS WITH ALMA 

Christopher Bradshaw 

Department of Astronomy, Yale University, New Haven, CT 06511 


Stella S. R. Offner 

Department of Astronomy, Yale University, New Haven, CT 06511; Department of Astronomy, University of Massachusetts, Amherst, 

MA 01003. 


Hector G. Arce 

Department of Astronomy, Yale University, New Haven, CT 06511 
Draft version February 2, 2015 

ABSTRACT 

Protostellar outflows provide a means to probe the accretion process of forming stars and their ability 
to inject energy into their surroundings. However, conclusions based on outflow observations depend 
upon the degree of accuracy with which their properties can be estimated. We examine the quality 
of ALMA observations of protostellar outflows by producing synthetic CO(1-0) and CO(1-0) ob¬ 
servations of numerical simulations. We use various ALMA configurations, observational parameters, 
and outflow inclinations to assess how accurately different assumptions and setups can recover un¬ 
derlying properties. We find that more compact arrays and longer observing times can improve the 
mass and momentum recovery by a factor of two. During the first ^0.3 Myr of evolution, ^^CO(l- 
0) is optically thick, even for velocities |u| > 1km s“^, and outflow mass is severely underestimated 
without an optical depth correction. Likewise, CO(1-0) is optically thick during the first 0.1 Myr. 
However, underestimation due to shorter observing time, missing flux, and optical depth are partially 
offset by the assumption of LTE and higher excitation temperatures. Overall, we expect that full 
ALMA CO(1-0) observations of protostellar sources within 500 pc with observing times > 1 hrs and 
assumed excitation temperatures of T < 20K will reliably measure mass and line-of-sight momentum 
to within 20%. 

Subject headings: stars: formation, stars:low-mass, stars:winds, outflows, ISM: jets and outflows, 
interferometry 


1. INTRODUCTION 


Young protostars are obser ved to launch energetic colli - 


mated bipolar mass outflows (Lada 1985; Bachiller 1996) 
These protostellar outflows play a fundamental role in 
the star formation process on a variety of scales. On 
sub-pc scales they entrain and unbind core gas, thus set- 


(iMatzner & McKee 

19991 lAlves et al. 
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Hosokawa||20E| OiJ 

mer Arce 2014p 

. Interaction be- 


stellar accretion and, ult imately, terminate it (Arce et al. 
2007 Frank et al. |2014). On sub-pc up to cloud scales, 
outflows inject substantial energy into their surround¬ 
ings, potentially providing a means of sus taining cloud 
turbulence over multiple dynamical times (Nakamura fc 
Li|2007| ICarroll et al.|2009||Wang et al.|2010||Arce et al 


20101 INakamura et al. 120111 Hansen et al. 2012 |Naka- 

- Lij2Q14^. -- -- 


I'he origin of outflows is attributed to the presence of 
magnetic fields, and a variety of different models have 
been proposed to explain the launching mechanism. Of 
these, the “disk-wind” model ( Blandford fc Payne||1982 


Pelletier & Pudritz 1992), in which the gas is centritu- 
gally accelerated from the accreti on disk surface, and the 
“X-wind” model (Shu et al. 1988), in which gas is acceler¬ 
ated along tightly wound field lines, are most commonly 


invoked to explain observed outflow signatures. However, 
investigating the launching mechanism is challenging be¬ 
cause launching occurs on scales of a few stellar radii and 
during times when the protostar is heavily extincted by 
its natal gas. Consequently, separating outflow gas from 
accreting core gas, discriminating between models, and 
determining fundamental outflow properties are nontriv¬ 
ial. 

Single-dish molecular line observations have been suc¬ 
cessful in mapping the extent of outf lows and their kine- 
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ever, outflow gas with velocities comparable to the cloud 
turbulent velocity can only be ex tracted with additional 
assumptions and modeling (e.g., 


Bontemps et al. 1996 


Arce fc Goodman|2QQl[[Maury et aI.|2UU9[ Dunhani et al. 


2014), which are difficult to apply to confused, clu^ 

terea star forming environments (e.g ^^ mtis et al |2Q1Q 


Graves et al. 2010 Plunkett et al. 2013|). Meanwhile' 
small scale structure and the inner launching region, 
which can provide clues about entrainment and vari¬ 
ability, are unresolved. Most interferometer observations 
have provided outflow maps with resolution of 10^ to 10^ 
AU scales (e.g.. 


Arce & Sargent 


Gueth & Guilloteau|1999 

2006||Hirano et al.|2{jl0| ), 


Lee et al. 2002 
feut a few recent 
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interferometer maps, including the latest maps obtained 
with the Atacama Large Millimeter/submilllimeter An¬ 
tenna ( ALMA) can probe down to scales of tens to about 
100 AU (iCodella et al.|2007l|Lee et al.|2009l|Cabrit et al 

-- ftflLeeet aL|20Ml|Oya € 


2012||Codella et al.|2014||Lee et al.|2014||Oya et al.|2014p 

interferometry is sometimes able to produce large high- 

resolution maps ([Plunkett et al. | |2013 Fernandez-Lopez 


et al. 2014 Storm et al. ' 2014[ ), but it can also resolve 


out larger scale structure if short-spacing or single-dish 
observations are not combined with the interferometer 
data. Consequently, it is sometimes difficult to assemble 
a complete and multi-scale picture of outflow properties 
with these observations. Uncertain optical depth, varia¬ 
tion in excitation and abundance, and projection effects 
further complicate the extraction of physical information. 

In the absence of either multiple emission lines or trac¬ 
ers, gas above a few km s~ ^ is often assumed to be opti¬ 
cally thin and in LTE (e.g. Cabrit fc Bertout|199Q | Lad a 
& Fich 1996 Bourke et al. 1997| Dunham et al.||2U14)' 


Analytic models provide a means to evaluate these as¬ 
sumptions, however most models assume outflow sym¬ 
metry, smooth velocity and densi ty distributions, and 


a simple temperature pro f ile (e.g., [Canto fc Raga||1991 
Masson fc Chernmll19931 Raga & Cabrit] |1993| |CantQ| 
et al.||2QQQ| ~ Lee et al.||2QQl Carolan et al.||2QQ8 ). In con¬ 


trast, three-dmierisioiialTnyarodyriarh^^ sup¬ 

ply detailed predictions for physical quantities related 


et al.j 

|2012t 

[Machida & Hosokawa[ 

|2013t 

[Oftner ^ Arce 


rigorously. Thus, the most promising avenue for con¬ 
structing a complete picture of outflows lies at the in¬ 
tersection of numerical modeling and observations. By 
performing synthetic observations to model molecular 
and atomic lines, continuum, and observational effects, 
simulations can be mapped into the observational do¬ 
main wher e they can be compared directly to observa - 
tions (e.g., Offner et al.||2Qll 2012a Mairs et al. 2014). 


Such direct comparisons are important for assessing the 
“reality” of the simulations, to interpret ob servational 


data and to assess observational uncertainties (Goodman 
2011). In addition to observational instrument limita¬ 


tions, chemistry and radiative transfer introduce addi¬ 
tional uncertaint ies that are difficult to q uantify without 
realistic models (Beaumont et al. 2013). Synthetic ob¬ 
servations have previously been performed i n the context 
of un derstanding outflow opening angles (lOffner et "ah 


2011), observed morphology ( Peters et ah I |2014|), and 


impact on spectral energy distributions (jottner et al. 

'20i2bD. 


I'he immanent completion of ALMA provides further 
motivation for predictive synthetic observations. Al¬ 
though ALMA will have unprecedented sensitivity and 
resolution compared to existing instruments, by nature 
interferometry resolves out large-scale structure and dif¬ 
ferent configurations will be sensitive to different scales. 
Atmospheric noise and total observing time may also ef¬ 
fect the fidelity of th e data. Prev i ous syn thetic obser¬ 
vations performed by jOffner et al.| (|2012a|) suggest that 
the superior resolution of lull ALMA and the Atacama 
Compact Array (ACA) will be able to resolve cor e struc¬ 


ture and frag mentation prior to binary formation. Peters 
et al. (2014) predict that ALMA will be able to resolve 
complex outflow velocity structure and helical structure 


in molecular emission. However, neither of these stud¬ 
ies was sufficiently comprehensive to assess the fidelity 
of ALMA observations under different assumptions and 
conditions. 

In this paper we seek to quantify the accuracy of 
various ALMA configurations in recovering fundamental 
molecular outflow properties such as mass, line-of-sight 
momentum, and energy. We use the CASA software pack¬ 
age to synthetically observe protostellar outflows in the 
radiat ion-hydrodynamic calculations of [Offner fc Arce| 
(2014). These simulations supply realistic distributions 
for gas densities, temperatures, and velocities, which 
more closely represent the complexity of observed out¬ 
flows than analytic models. By modeling the emission at 
different times, inclinations, molecular lines, and observ¬ 
ing configurations we evaluate how well physical quanti¬ 
ties can be measured in the star formation process. In 
section ^we describe our methods for modeling and ob¬ 
serving outflows. In section ^ we evaluate the effects 
of different observational parameters on bulk quantities. 
We discuss results and summarize conclusions in Q 

2. METHODS 

2.1. Hydrodynamic Simulation 

In order to assess the recovery of information us¬ 
ing ALMA, we post-process a self-gravitating radiation- 
hydrodynamic simulation in which we have complete 
three-dimensional temperature, density and velocity in¬ 
formation. The simulation we focus on here, thO.lfwO.3, 


was previously presented in (Offner & Arce 2014 here¬ 
after OA14). This calculation is not intended to model 
a specific source but rather to represent a “generic” 
molecular outflow originating from a low-mass protostel¬ 
lar source. The outflow is assumed to be driven by an 
X-wind with a collimation angle that results in a high- 
momentum region along the poles (which mimics a jet) 
and a lower momentum wider angle component. Since 
the simulation follows the evolution over 0.5 Myr it could 
in principle be used to model outflows from young pro¬ 
tostars (e.g., L1448, HH 212, HH 211) as well as slightly 
more evolved protostars (e.g., HH46/47, L1551-IRS5, 
L1228). We briefly describe the simulation properties 
below and refer the reader to OA14 for additional de¬ 
tails. 

The calculation was performed with the ORION adap ¬ 
tive mesh ref inement (AMR) code Truelove et al. (1998); 
Klein (1999). The simulation follows the collapse of 
an isolated, turbulent low-mass core. It begins with 
an initially uniform, cold AMq sphere of radius Rc = 


2 X 10^'cm, density pc = 2x 10“^^ g cm“'^ and tempera¬ 
ture Tc = 10 K. This core is embedded in a warm, diffuse 
gas with p = pc/100 and T = lOOTc = 1000 K. The dense 
gas is initialized with a grid of random velocity perturba¬ 
tions such that the initial rms velocity dispersion is 0.5 
km s“^. 

Additional levels of adaptive mesh refinement (AMR) 
are inserted as the core collapses under the influence of 
gravity. The core itself is resolved with a minimum cell 
size of Amin — 0.001 pc, where the maximum level of re¬ 
finement has Amin — 26 AU. Once the central region ex¬ 
ceeds the m aximum grid resolutio n (pmax — 6.5 x 10“^^ g 
cm“^, (e.g., Truelove et al.||l998)), a “star” forms. This 
star, which is represented by a Lagrangian sink parti- 
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cle, accr etes, radiates and launches a collimated bipolar 
outflow (iKrnmholz e t al.||2QQ4| |Offner et al.||2009| | Cun¬ 
ningham et al. 1120111. The rate of mass loss due to the 


ningham et al. pOTTlF 
outflow is set to a fixed fraction of the instantaneous ac¬ 
cretion rate: 171,^ = where = 0.3 is the outflow 

launc hing rate given by the X-wind model (Shu et al 


1988). The distribution of outflow momentum is param¬ 


eterized by a fixed collimation angle, Oq = 0 . 1 , which 
is empirically d etermined to be similar to that of ob 
served out flows (Matzner & McKee 1999 Cunningham 
et al. ]MTt . Although Uq is constant m time, the outflow 
injection into the AMR grid occurs on such small scales 
that the outflow properties such as the opening angle and 
morphology evolve hydr o dynamically. Previo us work by 
Offner et al. (2011) and Offner & Arce (2014) show that 
these properties agree well with observed low-mass out¬ 
flows. In particular, the outflow morphologies character¬ 
istic of CO e mission appear s imilar to those of observed 
outflows (e.g., Lee et al.|2002)and the simulations repro- 
du ce the outflow opening an^e dependence on time seen 
by Arce & Sargent (2006). 

In the simulation the protostar forms around t 0.17 
Myr, with the outflow launching beginning shortly there¬ 
after (when the protostellar mass reaches 0.05 Mq). For 
the first few kyr the outflow is confined to the core. It 
succeeds in breaking out around t ^0.22 Myr, at which 
point the gas exhibits well-defined “v-shaped” cavity 
walls. After breakout, the outflow entrains and ejects a 
significant fraction of the core, sweeping up progressively 
wider solid angles of gas. By t = 0.35 Myr, much of the 
core has been dispersed and the outflow structure is no 
longer clearly visible in the dense gas. Most of the resid¬ 
ual gas is either cold and clumpy or warm and diffuse. 
The protostar continues accreting from its local accretion 
disk, but infall has significantly diminished. The outflow 
continues but at two orders of magnitude below the ini¬ 
tial launching rate and little gas is entrained. As a result, 
the outflowing gas does not exhibit a clear signature in 
the cold (molecular) gas. 

By the end of the calculation (t = 0.5 Myr), the sim¬ 
ulation contains a single star with a mass of ^1.45 Mq. 
Most of the core mass has been ejected from the domain, 
and the remaining gas has a rms mass-weighted velocity 
dispersion of r\j 1 km s 

Since OA14 found that the final stellar mass and star 
formation efficiency did not depend strongly on Oq and 
we analyze only a single calculation. However, we 
note that different initial core masses, rotations, and 
magnetic field strengths might produce qu alitatively dif¬ 
ferent results (Machida & Hosokawa||2013). 

In the remainder of the paper, we focus on the colder 
molecular component of the simulated outflow, as would 
be traced by low-level rotational transitions of CO, rather 
than on the warmer jet component, which is predomi¬ 
nantly atomic or ionized gas here. The highest velocity 
gas (the jet) in the simulation reaches about 75 km s“^, 
which is similar to the deprojecte d flow speeds estimated 
for low-mas s sources like HH 212 (|Lee et al. |2008t |Cabrit 
et al.||2Ql^ . In contrast, the cold, molecular cornponent 
has much lower average velocities of a few km s“^. 


2 . 2 . Molecular Line Modeling 


We use the non-local thermodynamical equilibrium ra¬ 
diative transfer code radmc-3eH to compute the line 
emission in ^^CO(l-O) and ^^CO( 1 - 0 ). We ad opt the 
Large Velo city Gradient (LVG) approximation (Shetty 
et al. IMT] ), which solves for the rotational level popula- 
tions by solving the equations for local radiative statisti¬ 
cal equilibrium. RADMC-3D requires 3D input gas densi¬ 
ties, velocities and temperatures, which are produced as 
outputs by the hydrodynamic simulation. We perform 
the radiative transfer on a uniform 256^ grid, where we 
interpolate all the AMR data to the second refinement 
level (Ax = O.OOlpc). We include turbulent line broad¬ 
ening on scales at and below the grid resolution by adding 
a constant micro-turbulence of 0.05 km s“^. 

To obtain the CO abundances from the total gas den¬ 
sity, we assume that molecular Hydrogen dominates in 
all gas cooler than 1,000 K, where = p/{2.Smp). We 
adopt constant CO abundances of [^^CO/H 2 ]= 8.6 x 10 “^ 
(iFrerking et a l. 1982) and [^^CO/^^CO]=62 (jLanger & 
Fenzias 1993 )ior gas cooler than 900 K; otherwise the 
CO abundance is set to zero. Thus, we adopt a sim¬ 
ple prescription for the C O abundance and do not follow 
chem ical networks (e.g. Viti et al. 2004 Offner et al. 
2014). Consequently, the CO line smission only origi- 


nates from relatively cold gas in the dense core and gas 
entrained by the outflow; the warm, low-density ambi¬ 
ent material and the hot, outflow gas, which is ionized 
by construction, do not em it. We adopt the mo lecular 
collisional coefficients from Schoier et al. (2005). The 
output spectral cube has a velocity resolution of 0.08 km 
s“^and velocity range jxl < 10 km s“^. Little CO emis¬ 
sion originates from gas moving faster than 10 km s“^, 
because the simulated gas tends to be atomic or ionized. 
As cold material is swept up by the outflow, the net 
gas velocity diminishes. Observed outflows often have 
more emission at higher velocities than we find here (e.g 


Bachiller et al.||1991[ 1995] Cernicharo fc Reipurth||1996 


Cueth &: CuillQteau|1999| van der Marel et al.|2013[]Durr 


ham et al.||2014p. The difference could arise because the 

adopted launching model underestimates the amount of 
high-velocity gas or because the entrainment of cold gas 
is even more efficient than predicted by the simulations. 

2.3. Interferometry Modeling 

We convert the RADMC-3d outputs into skymaps 
in units of Janzky px“^ that we post-process us¬ 
ing the Common Astronomy Software Applications 
(CASAf] For our fiducial parameters, we adopt the 
distance, d = 450 pc, and sky position, (a, 5) = 
08^25^^415, -5R00'47" (J2000), of HH 46/47, an outflow 
recently observed with ALMA Cycle 0 (Arce et al.|2013 ). 
The simulation is not intended to be an exact model of 
HH46/47 but instead represent a typical young, isolated 
low-mass outflow of which HH46/47 is one example. 

We use the CASA task “simobserve” to produce a 
model observation for a given antenna configuration and 
set of observing conditions. We add thermal noise as¬ 
suming 0.5 mm of precipitable water vapor and a ground 
temperature of 269 K, which reflect good observing con¬ 
ditions. The synthetic maps are comprised of 27 point- 

^ http://www.ita.uni-heidelberg.de/ dullemond/software/radmc- 
3d/ 

^ http://casa.nrao.edu 
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ings, which cover an area of approximately 128” x 128”. 
We deconvolve and clean the synthetic observation using 
the task “simanalyze”. We vary the ALMA configura¬ 
tion, pointing time, cleaning parameters, and integration 
times as indicated in Table [TJ We find that the choice of 
cleaning parameters had little impact on the recovered 
mass (see the Appendix for further detail), and so we 
focus our analysis on the other parameters. 

2.4. Derivation of Mass, Momentum and Energy 

It is possible to estimate the mass from the emission 
a single CO line if several assumptions are applied 
Bourke et al.| (|1997|): (1) the gas is in local thermody- 
namic equilibrium (LTE), (2) the molecular line is op¬ 
tically thin, and (3) the level populations can be mod¬ 
eled with a single excitation temperature. One or more 
of these assumptions is routinely made when estirnating 
outflow properties in observations ( Dunham et al.|[2M4 , 
and references therein). When multiple lines or tracers 
are available a more rigorous analysis can be performed, 
however, constraining the accuracy of these approxima¬ 
tions without full information is challenging. By analyz¬ 
ing the simulations in which we have complete knowl¬ 
edge of the underlying mass and temperature distribu¬ 
tion, we can individually examine these approximations 
and quantitatively estimate the errors introduced. 

The simplest mass derivation adopts all three assump¬ 
tions. In this case, the outflow mass derived from a 
voxel with channel velocity v and position {a, (5) is given 
by: 

Mi2{v,a,S) = [R2/^^CO]ii^A{a,5)F{T,^)TBiv,a,5)Av 

( 1 ) 


where 


F(Tex) = 2.31 X 10 


14 


Tex + 0.92 


1 _ e-5.53/Te. j(2;j _ 


and 


J(T) = 


hu 


k{ekT — 1 ) 


( 2 ) 


( 3 ) 


Here, [H 2 /^^CO] is the abundance ratio of H 2 to ^^CO, 
l^rn = 2mp is the mean molecular mass of H 2 A is 
the area of emission (Ax^), Tex is the excitation tem¬ 
perature, Tbg is the background temperature, Tb the 
brightness temperature, v is the velocity, k is the Boltz¬ 
mann constant, u is the frequency of the line transi¬ 
tion (m2co(i-o)=113.3 GHz, z/i3co(i-o)=110-2 GHz), 
and h is Planck’s constant. Likewise, the mass es¬ 
timated from the ^^GO emission is Mis{v,a,6) = 
[^^GO/^^GO]Mi 3 (i;, a, S), where we adopt an isotopic ra¬ 
tio of [^^GO/^^GO] = 62. In observations, Tbg = 2.73 
K, however, we set Tbg = 0 K, since we do not include 
the cosmic microwave background in the modeling. In 
the fiducial case, we adopt an excitation temperature of 
Tex = 11 K, which reflects that the molecular outflow 
is dominat ed by entrained mat erial. This is the value 
adopted by Bourke et al. (1997), and we find that 11 K 


^ Note that typically, but not always, jUm takes into consid¬ 
eration the typical Helium and metal abundances in the cloud to 
obtain the total gas mass, in which case it has a value of 2.8mn for 
gas composed of 71% H, 27% He and 2% metals, see, e.g.,|Dunham| 
|et al.|2014| 


is in fact similar to the typical molecular gas temperature 
in the simulatio n (se e Table . We assess the accuracy 
of this value in §3.2| 

We derive the brightness temperature, T^, in each 
voxel from the emission computed by radmc-3d assum¬ 
ing the gas is a blackbody in the Rayleigh-Jeans limit: 


Tb = 




( 4 ) 


2v‘^k ’ 

where By is the intensity. The total mass is then 

rVmax r 

Mtot = / / M{v, a, 6)dvdad6, ( 5 ) 

Jvmin -'a,S 

where i;niin and i^max are either the expected velocity 
range of the outflow or the limits of the spectral cube 
(±10km s“^). 

To obtain the momentum and energy, we use the es¬ 
timated mass for each voxel as well as the line-of-sight 
velocity: 


^tot — 
^tot = 


/•r^max r 

/ / M{v, a, 6)vdvdad6, (6) 

dvmin da,d 

/ / -M{v, a, 6)v^dvdad5. (7) 

Jvmin ^ 


3. RESULTS 
3.1. Mass versus Time 

We observed the simulation at seven different times 
distributed from 0.2 Myr to 0.5 Myr, which spans the 
main accretion period of the forming protostar. In the 
fiducial case the outflow is inclined 20 — 30° with respect 
to the line of sight. We define the outflow as any emit¬ 
ting gas with velocity |i;| > 1km s~^ along the line of 
sight. This is consistent with observational determina¬ 
tions of the outflow mass from molecular emission which 


Hatchell et al. 20071 |Arce 

! et al. 20101 van der Marel 

et al. |2013t Dunham et al. 

2014p. This is the minimum 


velocity for which the outflow can be reliably separated 


from the embedding turbulent cloud (Arce & Goodman 
)02). In the simulation, the turbulent rms core velocity 
is ^ 0.5 km s~^ at the time of the initial outflow launch. 

Under this outflow definition, the derived mass initially 
increases as material is swept up and accelerated, reaches 
a peak around 0.25 Myr, and then declines when most 
of the core has been accreted. The decline is due to a 
combination of diminishing accretion, and hence outflow 
launching, and less material available for entrainment. 
We stress that this definition, which follows observational 
conventions, does not directly measure the full protestel- 
lar wind information, but rather the emission from the 
GO in the cloud entrained by the protostellar wind. 

Figure [l] shows the ^^GO outflow mass for velocities 
I'T’I > 1,1.5,2 km s~^. All exhibit similar trends of mass 
versus time. At early times the derived mass overesti¬ 
mates the underlying gas mass. This is mainly due to the 
assumption of LTE, which overestimated the line emis¬ 
sion. We discuss this further in the following section. 

Figure shows the fraction of the total mass esti¬ 
mated from ^^GO as a function of velocity. Approxi¬ 
mately 20% of the mass has a line-of-sight (los) velocity 
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TABLE 1 

CASA Run Parameters 


Run^ 

CO Isotope 

^los 

Conflguration 

hot 

tpointing 

-^clean 

-Ulean (^Jy) 

R1 

12 

30 

Compact 

7200 

10 

10000 

0.1 

R2 

12 

30 

Cycle 1.3 

7200 

10 

10000 

0.1 

R3 

12 

30 

Full 3 

3600 

10 

10000 

0.1 

R4 

12 

30 

Full 3 

7200 

10 

10000 

0.1 

R5 

12 

45 

Compact 

7200 

10 

10000 

0.1 

R6 

12 

45 

Cycle 1.1 

7200 

10 

10000 

0.1 

R7 

12 

45 

Full 3 

3600 

10 

10000 

0.1 

R8 

12 

45 

Full 3 

7200 

10 

10000 

0.1 

R9 

13 

45 

Compact 

7200 

10 

10000 

0.1 

RIO 

13 

45 

Cycle 1.1 

7200 

10 

10000 

0.1 

Rll 

13 

45 

Full 3 

3600 

10 

10000 

0.1 

R12 

13 

45 

Full 3 

7200 

10 

10000 

0.1 

R13 

13 

30 

Cycle 1.1 

7200 

30 

10000 

0.01 

R14 

13 

30 

Cycle 1.1 

7200 

30 

10000 

0.1 

R15 

13 

30 

Cycle 1.1 

7200 

30 

10000 

1.0 

R16 

13 

30 

Cycle 1.3 

7200 

10 

10000 

0.1 

R17 

13 

30 

Cycle 1.3 

7200 

30 

10000 

0.1 

R18 

13 

30 

Compact 

7200 

10 

10000 

0.1 

R19 

13 

30 

Compact 

7200 

30 

10000 

0.1 

R20 

13 

30 

Compact 

7200 

100 

10000 

0.1 

R21 

13 

30 

Full 3 

3600 

10 

1000 

0.1 

R22 

13 

30 

Full 3 

3600 

10 

10000 

0.1 

R23 

13 

30 

Full 3 

3600 

10 

20000 

0.1 

R24 

13 

30 

Full 3 

300 

10 

10000 

0.1 

R25 

13 

30 

Full 3 

600 

10 

10000 

0.1 

R26 

13 

30 

Full 3 

7200 

10 

10000 

0.1 


^ Run name, CO isotopologue, viewing angle with respect to the outflow axis, ALMA antenna conflguration, total integration time, pointing 
time, number of cleaning iterations, and the cleaning threshold. 



Fig. 1. — Observed (solid) and simulated (dotted) outflow mass 
versus time with various cutoff velocities. The observed mass was 
derived using the flducial full ALMA conflguration (R26 in Table 
p^ . The actual mass of the simulation along the same sightline for 
each cutoff is indicated by the dotted lines. 

above 1 km and 5% has a velocity above 2km s~^. 
This fraction is somewhat sensitive to the observing con¬ 
figuration. Configurations that resolve the large-scale 
structure well (e.g., R26) recover more of the emission at 
lower velocities. All observations miss the very low ve¬ 
locity gas, which has the highest optical depth, but most 
configurations detect similar amounts of higher velocity 
gas. 

3.2. Excitation Temperature 

Observational derivations of the outflow mass typically 
adopt a single excitation temperature. Under the as¬ 



Fig. 2.— Percentage of all mass in the domain with velocity 
greater than various minima for two synthetic observations and 
the simulation in ^^CO. Observations were made at t = 0.27 Myr 
with the R26 and R18 conflgurations. 


sumption of LTE, the excitation temperature is identical 
to the gas temperature. In equation we adopt a fidu¬ 
cial value of Tex =11 K. However, the gas has a range of 
underlying temperatures and is not guaranteed to be in 
LTE. Eigurej^ illustrates the impact of the assumed ex¬ 
citation temperature on the mass estimate. The higher 
the assumed excitation temperature, the higher the esti¬ 
mated mass, so errors in the assumed temperature have 
a large impact on the accuracy of the mass estimate. As 
shown in Table we derive mean and median values of 
^ 20 K and 10 — 11 K, respectively, for the overall gas 
temperature. The mean outflow temperature is signif¬ 
icantly higher than the median, which indicates that a 
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TABLE 2 

Average gas temperature 


hun(Myr)^ 

rtot(K) 


Median r(K) 

Median T|^|>i(K) 

0.20 

28.22 

69.56 

11.25 

10.30 

0.24 

18.68 

41.33 

10.06 

10.07 

0.27 

17.80 

28.03 

10.08 

10.15 

0.30 

17.08 

30.15 

10.08 

10.13 

0.35 

19.46 

68.07 

10.10 

11.70 

0.40 

19.62 

62.57 

10.22 

48.96 

0.50 

22.84 

69.90 

11.52 

73.86 


^ Simulation time measured from the start of the calculation and 
mass-weighted average gas temperature and median temperature 
defined over all velocities and over only high velocities. The average 
temperature is computed as T = 7 where the sum is only 

taken over gas with T < 300 K or gas with T < 300 K and hj > 
1km s“^. 



minority of warm cells are dominating the average. The 
very low median outflow temperature underscores the ex¬ 
tent to which the outflow entrains cold core gas. The me¬ 
dian and mean outflow temperature are only similar at 
the last two times when nearly all the core gas has been 
accreted or expelled, leaving little cold gas for entrain¬ 
ment. Based on our simulation comparison between the 
observed and actual mass, equation |l] gives a more accu¬ 
rate mass assessment when a value close to the median 
temperature is adopted and systematically overestimate 
the actual outflow mass for larger values. 

The origin of the disagreement is the assumption of 
LTE. Using the ratio of the level populations computed 
by the RADMC LVG radiative transfer calculation, we 
can compute the effective LVG excitation temperature: 


^ ^ ^^-hulikTi^vo) 


ni 


9i 


( 8 ) 


where ni are the upper and lower level populations 
{u = 1, I = 0), Qu, gi are the statistical weights of 
these levels {gilg<d = 3), is the transition frequency 
{ 1 / = m3Qo(i-o))5 and Tlvg is the effective excitation 
temperature. In LTE, this temperature is equal to the 
gas temperature. Eigure|^ shows that the effective exci¬ 
tation temperature is systematically lower than the ac¬ 
tual gas temperature, and thus the gas is not in LTE. 
Much of this difference is due to the outflow, while the 
slower (core) gas is closer to LTE. Thus, adopting an ex¬ 
citation temperature that is comparable to or larger than 
the gas temperature causes equation to over-estimate 
the outflow mass. As a secondary check, we repeat the 
radiative transfer calculation under the assumption that 
the level populations are in LTE. We find that outflow 
mass estimated from equation for a 7200s observation 
at 0.24 Myr is ^10% lower than when the radiative trans¬ 
fer is performed assuming LVG, and consequently, nearly 
identical to the actual outflow mass. 

In a comparison be tween simulated ou tflows and CO 
synthetic observations |Peters et al.|(|2014|) also found that 
estimated physical quantities exceeded the actual values, 
often by factors of 2-3. They do not explore the origin of 
this discrepancy, however, they also adopt an equation 
similar to which also assumes LTE. Consequently, it 
is likely that the underlying level pwulations in their 
calculations are likewise not in LTE. □ 


^ They adopt a much warmer excitation temperature than what 


Fig. 3.— Mass with |i;| > 1 versus time calculated from 
simulated emission for various values Tex- The actual simulated 
mass along the same sight line for |i;| > 1 is indicated by the dot- 
dashed line (purple triangles). 



-20 -10 0 10 20 
Ilvg“Ilte (K) 

Fig. 4.— Temperature difference between the effective excitation 
temperature, Tlvg (solid blue), and the actual gas temperature, 
Tl,te (dotted, green), at time 0.24 Myr for cells with temperatures 
less than 300 K (solid, blue) and cells with temperature less than 
300 K and line-of-sight velocities greater than 1km s“^ (dotted, 
green). N is the number of cells in each bin. 


3.3. Integration Time 

Eigurej^ illustrates the effect of different total integra¬ 
tion times on the estimated mass. As expected, longer 
observations were able to recover more emission. The dif¬ 
ference between a 5 minute and 2 hour observation could 
be up to a factor of 2 in mass. In both cases, the struc¬ 
ture of the outflow cavity is visible, but more of the larger 
scale structure is resolved out in the shorter observation. 

we adopt here (50 K). However, their simulated gas is initially two 
times warmer and the mean total gas temperature is comparable 
to or exceeds the assumed excitation temperature in two of their 
three cases. 
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At late times, the 5 minute observation overestimates the 
mass due to the low signal-to-noise, which creates false 
signal (see Figure]^) In contrast, the 1 and 2 hour inte¬ 
gration times have very good signal-to-noise and recover 
much of the detailed outflow structure. 

At the earliest output time, there was little difference 
in the mass estimated by the different integration times 
since the gas is very optically thick. Even at later times, 
the total mass is significantly underestimated, since the 
emission is optically thick in the central region for the 
low velocity gas at all times. Consequently, at best half 
the total mass can be recovered from CO(1-0) even 
with long observing times. Otherwise, the gas with 
\v\ > 1km s“^, which is essentially the outflow, is suf¬ 
ficiently optically thin in ^^CO that a 1 hour full ALMA 
observation is able to recover nearly all of the emission. 





Fig. 5. — Top: (a) Mass versus time derived from observed 

with various integration times. The solid lines show all gas while 
the dashed lines indicate the mass with |n| > 1km s“^. Bottom: 
(b) Mass versus time derived from observed with various 

integration times for gas with |n| > 1km s“^. The observations 
are performed with the full ALMA configuration and 10 second 
pointings (R23-R26). The actual simulation mass is indicated by 
the purple solid line with star symbols, and the actual simulation 
mass with \v\ > 1km s“^ is indicated by the purple dashed line 
with star symbols. 



Fig. 6 .— Top: (a) emission contour plot, colored by veloc¬ 

ity for t = 0.27Myr. The synthetic observation uses the full ALMA 
configuration 3, 300s integration time, and 10 second pointings 
(R24). Bottom: (b) emission contour plot, colored by veloc¬ 

ity for t = 0.27Myr. The synthetic observation uses the full ALMA 
configuration 3, 7200s integration time, and 10 second pointings 
(R26). Velocity is calculated as the mass weighted average at 
each position. The first contour is at —3km s“^ with a step of 
0.1km s-L 


3.4. ALMA Configuration 

The choice of ALMA configuration has a significant 
effect on the quality of the observation. Table lists the 
four ALMA configurations in three different observing 
cycles that we compare. 

An outflow spanning a few O.lpc at a distance of 450 pc 
is well-mapped by a beam of a several arcseconds. Figure 
[T shows the mass estimated from the Cycle 1.1 and 1.3 
configurations. Configuration 1.3 (R17) with its longer 
baselines and smaller beam was significantly less effective 
than the Cycle 1.1 (R14) observation. As illustrated by 
Figure although the smaller beam can recover more 
small scale detail it resolves out most of the larger scale 
emission, including much of the outflow cavity structure. 
Consequently, the signal-to-noise is significantly lower. 

As expected, increasing the number of antennas fur¬ 
ther increases the emission recovery. For example, de¬ 
spite having a smaller effective beam, full ALMA in con¬ 
figuration 3 produces a very similar map to the Cycle 1.1 
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TABLE 3 

ALMA Configurations 


Configuration^ 


-Vantenna 

dsep (ll^) 

dmax (m) 

^min (^) 

0 (”) 

Area(m^) 

Cycle 0 Compact 


16 

56.9 

126 

18.3 

5.23 

7240 

Cycle 1.1 


32 

71.4 

166 

15.1 

3.99 

14500 

Cycle 1.3 


32 

164 

443 

21.3 

1.49 

14500 

Full 3 


50 

93.2 

260 

15.0 

2.54 

22600 


^ Antenna configuration, total number of antennas, average dish separation, maximum baseline, minimum baseline, effective beam resolu¬ 
tion, and total area of the various ALMA configurations. To get these. We use the antenna configuration files provided with CASA. 






Fig. 7.— Top: (a) Mass versus time for ^^CO observations per¬ 
formed with the Cycle 1.1 (blue circles) and Cycle 1.3 (red crosses) 
configurations and an integration time of 7200s with 30s pointings. 
The actual mass is indicated by the green triangles. Bottom: (b) 
Mass versus time of synthetic ALMA observations. The observa¬ 
tions were performed with the Full Array configuration 3 and the 
Cycle 0 Compact configuration with an integration time of 7200s 
and 10s pointings. The actual mass is indicated by the green tri¬ 
angles. 


Fig. 8. — Top: (a) ^^CO emission contour plot, colored by veloc¬ 
ity for t = 0.27Myr. The synthetic observation is performed with 
the 1.3 ALMA configuration, over 7200s with 30 second pointings 
(R17). Bottom: (b) ^^CO emission contour plot, colored by ve¬ 
locity for t = 0.27Myr. The synthetic observation is performed 
with the 1.1 ALMA configuration over 7200s with 30 second point¬ 
ings (R15). Velocity is calculated as the mass weighted average 
at each position. The first contour is at —3km s“^ with a step of 
0.1km s“^ 


configuration map (compare Figure [sb with Figure]^ ). 
As illustrated by Figure [^, full ALMA performs signifi¬ 
cantly better with higher signal-to-noise than Cycle 0 for 
the same integration and pointing times, and a factor of 
^2 more mass is recovered. 

3.5. ^^CO(l-O) 

is often used to estimate outflow mass. While 
CO(1-0) is generally very optically thick in dense cores, 


it is less optically thick over the higher velocity outflow 
gas (e.g, |i;| > 2km s“^), and thus the assum ption of 
optical thinness and LTE is often ado pted (e.g. Lada & 
Fich||1996 van der Marel et ^[2013 ). This assumption 


IS more likely to be valid for low-resolution observations, 
in which the unresolved emission is spread over a larger 
area by the beam. However, at the high-resolution of 
ALMA we find that even the gas above a couple km s“^is 
optically thick in CO(1-0). Consequently, without an 







































Predictions for Observing Protostellar Outflows with ALMA 


9 


opacity correction the ^^CO severely underestimates the 
outflow mass. Figure shows the mass versus time for 
observations in ^^CO and ^^CO if we assume that both 
are optically thin. These results clearly demonstrate that 
^^CO is optically thick, even for gas with \v\ > 1. 



Time (Myr) 


Fig. 9. — Mass versus time for synthetic and obser¬ 

vations (R26, R24). The synthetic observations are performed with 
full ALMA configuration 3, 10 second pointings and 7200s total ob¬ 
servation time. The green triangles indicate the actual simulation 
mass. 


However, it is possible to use the optically thin 
emission to estimate the optical depth and then correct 
the mass estimate. Given the brightness temper¬ 

atures for each line, the relation between the bright¬ 
ness tempera tures and the optica l depth, ri 2 , can be 
expressed as (Dunham et al. 


1 - ^ [^^CO] 

T12 ?mbl3 [^^CO] 


( 9 ) 


riScoi 

where |i 2 QQj is the isotopic ratio. We can then correct 

the data by estimating the optical depth of each 

voxel and multiplying the emission in this velocity chan¬ 
nel by the appropriate factor. Figure pT| shows the esti¬ 
mated mass after this correction is applied. 

Figure shows the pre- and post-opacity corrected 
emission maps for ^^CO. Without the correction much 
of the outflow structure is missing simply because it has 
a low signal-to-noise. Applying the correction factor sys¬ 
tematically boosts the emission from the relatively dense 
outflow cavity walls, and the resulting map looks much 
more similar to that of map (e.g. bottom panel 

of Figure]^. Consequently, the is only useful for 

determining the velocity and structure of the outflow (at 
high resolution) when it is possible to measure the optical 
depth and apply a correction factor. The observa¬ 

tions, which in most cases do not require a correction, 
are significantly more accurate. 


3.6. Viewing Angle 

The accuracy of the mass recovery depends partially on 
the angle of the outflow axis along the line-of-sight. Here, 
we compare two different viewing angles. The fiducial 



Fig. 10.— Mass versus time for synthetic and ob¬ 
servation^ R2 6, R24), where the emission is corrected using 

equation The synthetic observations are performed with full 
ALMA configuration 3, 10 second pointings and 7200s total ob¬ 
servation time. The green triangles indicate the actual simulation 
mass. 


case is inclined 20 — 30° towards the line-of-sight and the 
second is inclined by 35 — 45° (the outflow axis is not 
fixed in the calculation and it varies by ^10°). 

While we expect slightly different results for the mass 
for different views. Figureshows that these particular 
sight-lines yield similar estimates. As long as the molec¬ 
ular emission is optically thin, as in the case of for 

\v\ >1 km s“^, the estimated mass should be similar 
for outflow inclinations from 15 — 75°. At extreme in¬ 
clinations the mass is far harder to estimate. Outflows 
viewed perpendicular to the axis (0°) are difficult to sep¬ 
arate from the ambient material because only a small 
component of the velocity is projected along the line-of- 
sight. Outflows viewed down the axis (90°) are easier 
to identify since the highest component of the velocity is 
along the viewing direction. However, at this inclination 
the two outflow lobes are projected on top of one another 
complicating the measurement of the gas morphology. 


3.7. Momentum and Energy 

The primary observational uncertainty in determining 
the line-of-sight momentum and energy relates to how 
well the mass can be recovered from the emission. In 
the case of an optically thin (or optical-depth corrected) 
tracer, we find that ALMA will be able to recover the 
mass and outflow structure quite well for reasonable ob¬ 
serving times and resolutions. However, estimating the 
total outflow momentum and energy is more fraught for 
several reasons. First, only one velocity component is 
directly observable, and the other two must be inferred 
using projection arguments and estimates of the incli¬ 
nation. Second, a significant fraction of the momentum 
may reside in lower velocity entrained material, which is 
both optically thick and difficult to distinguish from the 
surrounding dense core. If the source is sufficiently iso¬ 
lated, it is possibl e to model and include this material in 
the estimate (e.g. Arce fc Goodman 2001 Offner et al. 
2011 Dunham et al. | |2014| , however, accuracy relies on 
the assumed gas distribution of the parent core. Finally, 
outflow momentum and energy can reside in ionized or 
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TABLE 4 

Derived Physical Parameters^ 



^sim,tot 

■^sim, 1 v| >1 

^13CO,tot 


^l2co,tot 

^12 qo,|v|>1 

^12co,tot,corr 

^12co,|^;|>l,corr 

M (Mq) 

3.22 

0.39 

2.04 

0.49 

0.10 

0.051 

2.07 

0.50 

P (Mq km s 

1.57 

0.59 

1.38 

0.79 

0.12 

0.10 

1.39 

0.80 

E (10^2 erg) 

18.75 

13.45 

19.75 

16.47 

2.95 

2.84 

19.98 

16.67 


^ Mass, momentum, and energy for different mass definitions. All emission estimates are calculated with full ALMA configuration 3 (R26) 
with an integration time of 7200s and a pointing time of 10s for trun=0.27Myr. 


TABLE 5 
Derived Mass^ 


Run (Myr) 

Mtot 

^|u|>l 

Misco 

^13CO,|u|>l 

Mi2(JO 

^12 qo,|u|>1 

^12co,corr 

■^12cO, |u| >1,corr 

0.20 

4.22 

0.18 

1.70 

0.28 

0.08 

0.03 

1.71 

0.28 

0.24 

3.75 

0.47 

1.84 

0.53 

0.10 

0.06 

1.87 

0.54 

0.27 

3.26 

0.39 

2.04 

0.49 

0.10 

0.05 

2.06 

0.50 

0.30 

2.76 

0.29 

1.718 

0.40 

0.11 

0.04 

1.74 

0.40 

0.35 

1.89 

0.06 

1.01 

0.04 

0.09 

0.03 

1.10 

0.04 

0.40 

1.49 

0.04 

0.84 

0.02 

0.07 

0.01 

0.84 

0.02 

0.50 

1.31 

0.04 

0.74 

0.03 

0.11 

0.02 

0.74 

0.03 


^ Output time, total simulated outflow mass, simulated outflow mass for \vz\ > 1, total mass estimated from ^^CO, mass estimated from 
for \vz\ > 1, total mass estimated from with no opacity correction, mass estimated from for \vz\ > 1 with no opacity 

correction, total mass estimated from with opacity correction, mass estimated from for |i;^| > 1 with opacity correction. All 

emission estimates are calculated with full ALMA configuration 3 (R26) with an integration time of 7200s and a pointing time of 10s for 
trun=0.27Myr. All masses are in units of Mq. 


TABLE 6 

Derived Momentum^ 


Run (Myr) 1 

Psini,tot 

^sim,\vz 1 >1 

^13CO 

^13cO,|vz|>l 

^12CO 

^l2co,|vz|>l 

Pi2(JqCOTT 


0.2 

1.40 

0.26 

1.10 

0.60 

0.13 

0.11 

1.11 

0.60 

0.24 

1.85 

0.79 

1.44 

0.92 

0.15 

0.13 

1.46 

0.93 

0.27 

1.56 

0.59 

1.36 

0.79 

0.12 

0.10 

1.39 

0.80 

0.3 

1.33 

0.42 

1.11 

0.59 

0.10 

0.07 

1.12 

0.60 

0.35 

0.76 

0.11 

0.50 

0.10 

0.07 

0.04 

0.50 

0.10 

0.4 

0.52 

0.07 

0.35 

0.09 

0.04 

0.02 

0.35 

0.10 

0.5 

0.43 

0.07 

0.31 

0.10 

0.07 

0.03 

0.31 

0.10 


^ Output time, total simulated outflow momentum, simulated outflow momentum for \vz\ > 1, total momentum estimated from ^^CO, 
momentum estimated from for \vz\ > 1, total momentum estimated from with no opacity correction, momentum estimated 

from for \vz\ > 1 with no opacity correction, total momentum estimated from with opacity correction, momentum estimated 

from for \vz\ > 1 with opacity correction. All emission estimates are calculated with full ALMA configuration 3 (R26) with an 

integration time of 7200s and a pointing time of 10s for trun=0.27Myr. All momenta are in units of MQpms”^. 


TABLE 7 
Derived Energy^ 


Run (Myr) 1 

-R'sim,tot 

-‘^sim,\vz 1 >1 

Eisco 

^13CO|vz|>l 

Ei2qo 

^12qo|Vz|>1 

Ei2(JQ,COTT 

^i^coivz^l.corr 

0.2 

12.43 

6.98 

22.16 

19.58 

6.16 

6.06 

22.34 

19.73 

0.24 

23.12 

17.52 

22.91 

19.92 

3.19 

3.09 

23.19 

20.16 

0.27 

18.75 

13.44 

19.74 

16.47 

2.95 

2.84 

19.98 

16.66 

0.3 

15.65 

10.53 

13.40 

10.48 

1.47 

1.32 

13.56 

10.61 

0.35 

7.15 

3.85 

6.435 

4.31 

1.10 

0.96 

6.51 

4.36 

0.4 

4.20 

2.08 

6.67 

5.52 

0.48 

0.35 

6.74 

5.58 

0.5 

3.68 

1.99 

6.84 

5.89 

0.70 

0.51 

6.91 

5.96 


^ Output time, total simulated outflow energy, simulated outflow energy for \vz\ > 1, total energy estimated from ^^CO, energy estimated 
from for \vz\ > 1, total energy estimated from with no opacity correction, energy estimated from for \vz\ > 1 with 

no opacity correction, total energy estimated from with opacity correction, energy estimated from for \vz\ > 1 with opacity 

correction. Energy is in units of 10^^ erg. All emission estimates are calculated with full ALMA configuration 3 (R26) with an integration 
time of 7200s and a pointing time of 10s for trun=0.27Myr. 
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Fig. 11.— Top: (a) Uncorrected emission contour plot, 

colored by velocity for t = O.SMyr. Bottom: (b) Corrected (using 
equation emission contour plot, colored by velocity for 

t = O.SMyr. The observations use the full ALMA configuration 3, 
7200s integration time, and 10 second pointings (R4). Velocity is 
calculated as the mass weighted average at each position. The first 
contour is at —3km s“^ with a step of 0.1km s“^. 


very low-abundance warm gas, which is difficult to de¬ 
tect. In the simulations we analyze, the highest velocity 
material does not emit in CO, a situation which worsens 
at late times when there is little residual cold gas avail¬ 
able for entrainment. This gas contains a small amount 
of the outflow mass but a significant fraction of the out¬ 
flow momentum and energy. Here, we simply compare 
the line-of-sight momentum and energy, with the caveat 
that the total amount of momentum and energy account¬ 
ing for these factors is much higher. 

Tables and [^compare the simulation momentum and 
energy for a particular sightline with that derived from 
the observations. Tables and represent best-case 
ALMA observations with good mass recovery. The total 
momentum is underestimated by ^ 30% when all veloci¬ 
ties are considered, while the total energy is comparable 
or higher than the simulated value. This suggests that 
even applying a model for the core gas could underes¬ 
timate the outflow momentum if the mass is underesti¬ 
mated. The momentum and energy of the \v\ > 1 km s~^ 
gas is over-estimated due to the LTE approximation as 



Fig. 12 .— Mass versus time for synthetic observations for 

two views for inclinations differing by ~ 15 — 20°. The synthetic 
observations use the full ALMA configuration 3, 7200s total obser¬ 
vation time, and 10 second pointings. 



Time (Myr) 


Fig. 13.— Mass |n| > km s“^versus time for synthetic 
observations for two outfiow inclinations. The synthetic observa¬ 
tions are performed with the full ALMA configuration 3, 10 second 
pointings and 7200s total observation time. 


discussed in ^3.2 


Once corrected for opacity, observations give 

very similar momentum and energy estimates to the 
ones. Without the correction, observations 

significantly underestimate both by factors of 2-10. 


4. CONCLUSIONS 

We performed synthetic ALMA molecular line observa¬ 
tions of a numerical simulation of a forming protostar in 
order to evaluate the accuracy of inferred physical quan¬ 
tities. To produce the synthetic observations we adopt 
the distance of source HH46/47, which has recently been 
observed in CO by ALMA. Given it’s distance of 450 
pc and molecular extent of a few tenths of a parsec, it 
is a reasonable representative of isolated, local low-mass 
protostellar outflows. 

Following observational convention we adopted a con¬ 
stant excitation temperature of 11 K and assumed the 
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gas was in local thermodynamic equilibrium (LTE). We 
demonstrated that higher assumed excitation tempera¬ 
tures can easily increase the apparent mass by factors 
of 2-3. Inspection of the level populations of the out¬ 
flow gas (as defined by gas with line-of-sight velocity 
\v\ > 1km s~ ) indicate that the gas is not in LTE and 
that the effective excitation temperature is slightly lower 
than the actual gas temperature. This results in an over¬ 
estimate of the outflow mass, momentum and energy. 

By comparing with the actual simulated values, we 
found that the CO(1-0) line is an inaccurate tracer at 
high resolution since it is optically thick for all times 
and for velocities greater than 1km s“^. In contrast, 
the higher velocity gas emitting in CO(1-0) is optically 
thin for nearly all times. Consequently, the synthetic 
^^CO observations reliably recovered the gas mass for a 
variety of angular resolutions, viewing angles, and inte¬ 
gration times. Shorter pointing times resulted in slightly 
better results, but the difference was minor compared 
to the total integration time. The ^^CO mass estimates 
were significantly improved by applying an optical depth 
correction factor computed from the CO(1-0) line. 

Although full ALMA has the ability to achieve sub- 
arcsecond resolution, we found that significant flux was 
resolved out for effective beam sizes smaller than ^ 2 
arcseconds. Adding more antennas and using longer in¬ 
tegration times increased the estimated mass accuracy 
and structure recovery significantly. One to two hour in¬ 
tegration times recovered twice as much flux (and mass) 
as 5 and 15 minute observations. The best ALMA con¬ 
figuration we tested (full ALMA, configuration 3), was 
able to produce a detailed and high signal-to-noise map 


of the outflow with sufficient flux recovery that it over¬ 
estimated the outflow mass by ^ 20%. This was due to 
the assumption of LTE, which enhanced the conversion 
of flux to mass. 

The momentum and energy of the |'r’| > 1km s~^ gas 
were comparable to or greater than the simulated values 
for the best ALMA resolution. This is encouraging for 
observational assessments of the turbulent energy injec¬ 
tion due to outflows. However, corrections to account 
for outflow orientation and optical depth (for ^^CO) will 
be necessary to produce high accuracy estimates of the 
outflow impact. 

In summary, our analysis indicates that ALMA will al¬ 
low unprecedented study of molecular outflows with high 
resolution, accuracy and detail. However, additional un¬ 
certainties, including the molecular abundance, fraction 
of low-velocity material, and mass of the non-molecular 
component, will remain complicating factors even in the 
era of ALMA. 
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APPENDIX 




Time (Myr) Time (Myr) 

Fig. 14.— Left: (a) Mass versus time for synthetic observations with different numbers of cleaning cycles. The synthetic observations 

use the full ALMA configuration 3, with a 3600s integration time and 10s pointings (R21, R22, R23). Right: (b) Mass versus time for 
synthetic with various pointing times. The synthetic observations use the Cycle 0 compact ALMA configuration with a 7200s 

integration time (R18, R19, R20). 


There were several additional parameters that we explored but which turned out to have minimal impact. We discuss 
these parameters here for completeness. The parameter, Ydean, sets the limit for the number of cleaning iterations 
performed by “simanalyze”. We found that changes in the number of iterations had a small affect on the results. 
Eigurep^ shows the estimated mass for different values of Ydean- We adopted Ydean =10,000 as it seems to provide 
a balance between good results and efficient use of time and processor cycles. 
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The parameter Fdean sets the threshold for early termination of the cleaning. Our calculations never satisfied the 
early termination for any of 0.01 — O.lmJy, so we adopted a default value of 0.1 mJy. 

We also considered the effect of different pointing times. The CASA documentation suggests that the pointing time 
could have a significant effect on the observation, with shorter pointings giving better results. However, we found 
that the difference between 10s pointings and 100s pointings was marginal. Figure [l^ shows that observations with 
shorter pointings perform slightly better, however the integrated emission maps appear nearly identical. We use the 
CASA recommended 10s pointing time as our default value, but we suggest that the choice of pointing time may have 
little impact on actual observations. 
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